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We have considered the processes which lead to elastic scattering between two 
far ultraviolet or X~ray photons while they propagate inside a solid, modeled as a 



' simple electron gas. The new ingredient, with respect to the standard theory of 

' photon-photon scattering in vacuum, is the presence of low-energy, nonrelativistic 

' electron-hole excitations. Owing to the existence of two-photon vertices, the seat- 
ed ■ 



tering processes in the metal are predominantly of second order, as opposed to fourth 
order for the vacuum case. The main processes in second order are dominated by 



' exchange of virtual plasmons between the two photons. For two photons of similar 

2 ' energy HQ, this gives rise to a cross section rising like il^ up to maximum of around 

lO"'^'^ cni'^, and then decreasing like $7"^. The maximal cross section is found for the 

' photon wavevector k kp, the Fermi surface size, which typically means a photon 

\ 

energy in the keV range. Possible experiments aimed at checking the existence of 
these rare but seemingly measurable elastic photon-photon scattering processes are 
discussed, using in particular intense synchrotron sources. 
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I. INTRODUCTION 



Photon-photon scattering processes are well known in vacuum. The lowest order per- 
turbation theory diagrams are shown on Fig.|l]. The corresponding scattering amplitude and 
cross section, first calculated by Karplus and Neumann were later thoroughly investi- 
gated by de Tollis and are extensively reviewed in Lifshitz and Pitaevskii's book |Q. The 
cross section of the process in vacuum turns out to be very small, of the order of magnitude 
of lO^'^^cm^ at the threshold for the production of real — pairs in the intermediate 
state, which corresponds to an incident photon energy of hfl = mc^ ~ 0.5 MeV. For en- 
ergies hVt -C mc^, the cross section falls as a ~ {^^Y- It is thus extremely small in the 
eV range, where powerful laser beams are available, or even in the keV range, where strong 
synchrotron sources are now feasible. The physical reason is clear: if the energy of photons 
lies too far below the "energy gap" 2mc^, the virtual intermediate e~ — pairs are hard 
to excite, and this greatly cuts the probability of the process. In the end, this makes the 
process very difficult to pursue experimentally - albeit not impossible |jTl|, ||12 . 



One is thus led to the obvious idea that, if instead of vacuum, the photons were embedded 
in a material medium, where pair excitations could be created with a smaller or zero energy 
gap, then the probability of photon-photon scattering in such a system might be enhanced. 

Condensed matter physics provides us with countless examples of systems with this prop- 
erty: electron-hole pairs replace electron-positron pairs, and this reduces excitation energies 
from MeV to the eV or meV range. For example, in a semiconductor the intermediate state 
involves pairs excited across the semiconducting band gap. 

Third order nonlinear optical susceptibilities, which in principle involve the same di- 
agrams of Fig.0, have for example been well characterized, both experimentally 0, and 
theoretically 0, [|], for the quasi-one-dimensional prototypical semiconductor polyacety- 
lene. 

The focus of this paper will be however on the elastic scattering cross section among very 
high energy photons, which makes the problem radically different from those studied so far. 
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To make a start, we would like to study photon-photon scattering inside the simplest 
and most typical condensed system, for example that of a simple metal. Of course, even 
simple metals have in real life a number of complications. The periodic potential brings 
about a multi-band problem. Due to the multiplicity of bands, interband transitions may 
occur, etc. However we shall restrict in this paper to the contribution of intraband transitions, 
right at the Fermi surface, and these processes can be modelled with a single parabolic band. 
Moreover, in the keV photon energy range (which will turn out to be the most interesting 
case) virtually all systems, including insulators, can often be treated as free electron metals 
to a very good approximation. 

A second complication is that in a metal, or more generally in a solid or liquid, photons 
are heavily absorbed. Their energy decays into electron-hole pairs, and finally into heat. Of 
course, this fact may represent a practical obstacle at getting strong photon fields inside an 
actual material. Yet, for example, one can think of using a very intense evanescent photon 
field at a surface in a reflection experiment, or use intense, far ultraviolet and X-ray beams 
from a synchrotron radiation source, tuned so that the the absorption cofficient is small 
enough to be neglected. 

A third relevant issue might be the difficulty of separating elastic from inelastic photon- 
photon processes, when the process takes place inside a metal. Since there is no gap, a 
shower of soft electron-hole pairs is likely to accompany the photon-photon collision, and 
it could be difficult in practice to separate this inelastic part from the elastic process we are 
after. 

Postponing further discussion of some of these, and other related issues to a later section 
of this paper, we will now proceed with the actual calculation of the elastic photon-photon 
cross section in an electron gas, so as to come up with predictions for orders of magnitude 
and angular dependences. This calculation is not simply a generalization of the classic 
calculation of from vacuum to a uniform electron gas. The nonrelativistic nature of the 
electron gas allows for new processes, which do not exist in vacuum. There are new energy 
scales, the Fermi energy Ep and the plasma energy hup, which are absent in vacuum. 
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Anticipating our later results, we will end up with numerical values of the cross section, 
which are again exceedingly small - we estimate for a potassium target a peak value of 
10"'^^ cm^, somewhat smaller than that in vacuum. However, the peak is now in the keV 
photon frequency range. Since ultraviolet and soft X-ray photon sources are so much more 
intense than 7-ray sources, the possibility that our calculations might not remain academic 
seems real. Moreover, our estimated photon-photon cross section in the ultraviolet is ~ 
10~^^ cm^, which might be measurable using medium power lasers. At this stage, it is not 
clear to us whether such measurements are going to be feasible, or worth doing, or whether 
the present calculation will remain an academic exercise. We simply would like to point 
out that the new photon-photon processes exist, and that they are in principle finite and 
measurable. 

As announced, we shall take the uniform, single-band non-relativistic electron gas (EG) 
as our model for the metal. The model is completely specified by the electron mass m and 
the Fermi energy Ep. In section 2 we identify the simplest possible photon-photon scatter- 
ing processes, which now appear in the second, third and fourth order of the perturbation 
theory (as opposed to just fourth order in vacuum). In later sections, we choose two of these 
processes and calculate them out in detail. First, in section 3 we deal with the 2nd order 
processes, which are the largest, and consist essentially of exchange of plasmons between 
the photons. These processes do not exist in vacuum. In section 4 we then investigate the 
4th order process, formally the same square diagram of Fig.[^. Now, however, positrons are 
replaced by holes in the Fermi sea, which modifies both the form and the magnitude of the 
resulting scattering amplitude. The resulting fourth-order cross section, not surprisingly, 
turns out in the end to be many orders of magnitude smaller than that of the 2nd order pro- 
cesses, even at very high frequencies (where they could in principle have come closer). We 
therefore identify plasmon exchange as the dominant mechanism for the photon-photon scat- 
tering in the EG. In section 5 we summarize our results, and discuss briefly the possibilities 
of experimentally observing the proposed photon-photon elastic scattering processes. 
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II. POSSIBLE SCATTERING PROCESSES 



Before moving on to the actual calculation, we wish to stress some points, which make 
the situation of photons inside the EG different from that in vacuum. First of all, a photon 
in an electron gas is a well-propagating mode only at frequencies much larger than the 
plasma frequency cUp. Approaching that frequency from above, photons become dressed, 
absorption increases, until below cUp photons do not propagate anymore, and get absorbed 
within a wavelength or so. In an actual experimental situation, a beam from an external 
source penetrates into the metal, and the scattered, reflected or transmitted radiation, which 
leaves the metal, is detected. In these conditions, what actually scatters inside the metal, 
are the dressed photons, as they enter through the surface. To predict the outcome of 
such experiment rigorously, one would have to perform the calculation including photon 
dressing, plus the effect of the surfaces. Both represent nontrivial complications, which we 
shall not endeavor to consider at this stage. Therefore, we simply assume here free photons 
propagating in an infinite bulk electron gas. We shall restrict our treatment to frequencies 
n higher than Up 



and completely neglect the effects of photon dressing, and surface effects. This will allow us 
to use essentially the same formalism as in the case of vacuum. 

Since we want to treat the electrons as non-relativistic, we shall also restrict the energy 
hD, of our photons from above 



We shall also assume, trivially, that the Fermi energy of the electron gas (imagining typical 
metal density) satisfies 



(1) 



(2) 



(3) 



obviously satisfied in practice. 



5 



The starting point is the non-relativistic hamiltonian for electrons interacting with elec- 
tromagnetic field in the Coulomb gauge 



+ lJ J d^Xid^X2^^(xi)¥(x2)U(xi - X2)^(X2)^(XI) , (4) 

the electron charge being — e. Here ^'^(j;) creates an electron at point x (spin indices are 

J. — * — * 

omitted for simplicity) and al creates a photon of wavevector k and polarization A. The A 
operator is 



2\ 1/2 



where V is the normalization volume. The total hamiltonian can be split into the following 
terms 

H — Ho + Hcoul + Hrad + Hn + Hi2 , (6) 



where 



and 



Ho= f d^x'if\x)^f^{x) (7) 

J Zi lib 



Hrad = E ^^k4\^kX (8) 
kX 



are the free fields, 



Hcoul = \j j (i'x^d^X2^\x{)-^\x2)U{x^ - X2)*(X2)*(X1) (9) 



is the Coulomb interaction between the electrons and 



Hii^ — I d^x-^Ux)A(x).p<lf(x) (10) 
mc J 



and 

o2 



Hi2 = 2^ J d^x¥{x)^{x)A{x).A{x) (11) 
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are the interaction terms linear and quadratic in A{x), respectively. 

Now to have a non-zero S-matrix element Sfi between the initial photon state containing 
2 photons ki, k2 and the final state with two photons ^i, we must have in the S-operator 
S{—oo, +00) a term with at least four A operators 

Sfi = I {%h I a^'^S{-oo, +oo)at^ at^ \ %h) I $fs) , (12) 

I $p/j) and I being the photon vacuum and the ground state of the Fermi sea of electrons, 
respectively. The photon polarization index A has been absorbed into k. Elasticity of the 
process is implied by the Fermi sea being left in its ground state at the end of the process. 

For photon-photon scattering in vacuum, described by the relativistic hamiltonian, where 
the interaction term is linear in A, the lowest order of perturbation theory giving nonzero 
contribution is the 4th p|, 0]. The corresponding processes are the three diagrams on the 
Fig.|I|, differing from each other by the assignment of the external legs to the vertices of the 
square, plus three other diagrams which differ from the former ones just by the orientation 
of the internal fermionic loop. 

In our case, however, there are more possibilities, because the interaction term Hi2 is 
quadratic in A. The simplest scattering processes now appear in second, third and fourth 
order of perturbation theory, and are characterized by the respective scattering amplitudes 

r+00 



r+co 

J_ dt2 {<I>FS I {%h I a^- flfc- T[Ha{h)H,2{t2)] a\a\^ \ %h) I ^fs) , (13) 
f+00 

/ dts {^FS I (V I T[Hii{h)Ha{t2)H,2{t3)] al al \ V) I ^fs) , (14) 

J-00 12 Kl K2 



00 



f+00 

dh I (V I %X T[Ha{h)Ha{t2)Ha{h)Ha{U)] alal^ \ %h) \ , (15) 



where all operators are now in the interaction picture. In the following sections we shall 

?}^) and 
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investigate two of the above processes, Sf^ and in actual detail, and evaluate their 



scattering amplitudes. 

Since we shall finally be interested in the cross sections of the scattering processes, we 
recall that the usual scattering amplitudes Mfi satisfy the relations 

% = 6f, + t{27f)H{Pf - P,)Tf, , (16) 



Tfi = , . . . Mfi , (17) 

and that the differential scattering cross section ^ is related to Mfi through the equation 



da I \M 



fi 



|2 



(18) 



do QA'K'^c\ni + Vt2-c{ki + k2).nY' 

where fii,2 and ki^2 are frequencies and wavevectors of the incident photons and n is a unit 
vector in the direction of the solid-angle element do. 

III. SECOND ORDER SCATTERING PROCESSES: PLASMONS 

In this section we shall deal with the processes which result from the second order 
scattering amplitude (|l3l). In the end, these will be the dominant contribution to the 
scattering, due to the smallness of the fine structure constant a = ^ = Substituting 
for H.12 from ([lT|) and using the fact that the \1/ and A operators commute we get 

Sf = (y)'^(^)'/ d'x, J d'x2{^Ps\T[^Kx,Mxi)^Hx2Mx2)] \^FS) 

X {%h\ a,-a,^;T[(l(xi).l(xi))(l(x2).l(x2))]at 4^ | V) . (19) 

Using Wick's theorem, we find that there are six diagrams, three of which are shown on 
Fig.0. The three remaining diagrams differ from these just by interchange of the vertices 
and thus yield exactly the same contribution. Diagrams 1 and 2 of Fig.^ can be regarded as 
processes where the 4-momentum transfer between the photons is accomplished via exchange 
of an electron-hole pair. Diagram 3 describes instead the creation of an electron-hole pair 
by two-photon absorption and subsequent decay of the pair into two final photons. 
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Performing the usual algebra and introducing the electron Green's functions we get 
for Mj^-* from the bare diagrams on Fig.^ (including the factor of 2 from other 3 diagrams 
which just cancels the ^ factor) 

diagrams V / 

= (^)' (4vr/ic2)' E {ea-^.2)ie,3-eu)lfip,). (20) 

diagrams 

In the last expression we have introduced the usual Lindhard complex polarizability function 
n°(pj) of the free electron gas ^ , with pi for the i-th diagram determined by the 4- momentum 
conservation in the vertex. The unit polarization vectors of the photons incident with the 
vertex 1, resp. 2 of the i-th diagram have been denoted as en, and e^s, ej4, respectively. 

To proceed further, we notice an important difference between the interaction hamilto- 
nians Hn and Hi2. In a homogeneous medium, the electron- hole pair created by a photon 
through Hn is transverse and therefore cannot decay into a longitudinal Coulomb interac- 
tion. Hence, the process of Fig.^ has zero amplitude. If however, as in the present case 
(Fig.0), the electron-hole pairs are created by Hi2, this selection rule is absent, and there 
is nothing to prevent decay of the pair and a subsequent creation of another one due to 
Coulomb interaction. Therefore we must renormalize the bare diagrams of Fig.^ Within 
the random phase approximation (RPA) this is done by summing the whole "bubble" series 
of Fig.^, which we recognize as the familiar plasmon series PJ. The renormalized diagrams 
1 and 2 can now be interpreted as photons interacting with each other via an exchange of 
virtual plasmons. The diagram 3 then corresponds to absorption of two photons with cre- 
ation of an intermediate plasmon, plus subsequent decay of the plasmon into the two final 
photons. 

As is well known, exact summation of the plasmon series is equivalent to replacing the 
bare polarizability Il^{pi) by the screened polarizability 

Ti, 1 n°(p.) 
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where 

Uo{q) = ^ 

is the bare Coulomb interaction. Our final result for the scattering amplitude then reads 

diagrams 'jyt-ij \t-i-i 

Because the Lindhard function 11° (p) is itself a complicated function of the 4- momentum 
transfer p, the scattering cross section resulting from the last expression will also depend in 
an intricate way on the frequencies of the incident photons as well as on the geometry of 
the situation, which is determined by the wavevectors and polarizations of all the photons 
involved. In the rest of this section, we illustrate just a few particular cases. 

First of all, we notice that if the 4-momentum transfer is such that the denominator 
1 — f/o(pj)n°(pj) becomes zero, the scattering amplitude p^ ) diverges. This corresponds to 
an excitation of a real, instead of a virtual, plasmon. This unphysical divergence, however, 
is just a consequence of the random phase approximation, which we have used for simplicity. 
Summing a larger class of diagrams and dressing the plasmon would cause the pole to move 
away from the real axis, and thus remove the divergence. Nevertheless, we still expect the 
scattering cross section in a realistic metal to be considerably enhanced due to this "plasmon 
resonance" . 

Plasmon resonances can be hit as a function of scattering geometry. The result ( ^21) is 
valid for a generic geometry. For illustration, we have investigated simple planar scattering 
pictorially described in Fig.^. For simplicity, we further restrict the two incident photons 
to have equal energies and perpendicular wavevectors /ci ± A;2 (6' = 7r/2). The wavevectors 
of all four photons involved are therefore assumed to lie in the same plane, the polarization 
vectors being perpendicular to the latter. On Fig.^ we show the dependence of ^ on the 
angle /3 for h£lx = hVL2 = 10 eV. The density of the electron gas was taken to be that 
of metallic potassium, with kp = 0.73 x 10^ cm^^, corresponding to Ep = 2.03 eV. The 
corresponding plasma frequency is hup ~ 4 eV, well below hQi, hQ2, as required. One can 
see sharp peaks corresponding to the plasmon resonance. The peaks occur for angles (3 such 

10 



that the 4-momentum of the scattered photon, determined by conservation rules, exactly 
fits the requirement for excitation of a real plasmon in one of the diagrams. 

This calculation shows that, at least away from the peaks, the differential cross section 
for ultraviolet photons is exceedingly small, of the order of 10~^^ cm^. It is therefore of 
interest to ask what the dependence on the photon frequencies will be, and particularly 
whether the cross section could get larger in some other regime. In order to illustrate the 
frequency dependence, we consider the particular case of scattering of two incident photons 
with opposite wavevectors (Fig.^, 6 = n), assuming again all the photons involved to be 
polarized perpendicularly to the plane in which the wavevectors lie. In this case, in the 
diagram 3 of Fig.^ the 3-momentum transfer is zero, and this diagram does not contribute 
at all. In diagrams 1 and 2 the energy transfer is zero, and therefore no plasmon resonance 
is possible. It is easy now to find the limiting forms of ^ for cases k <^ kp (low frequencies) 
and k^ kp (high frequencies). 

For the case k <^ kp, we expand n°(pj)/[l — f/o(pj)n°(pj)] in (p2D in powers of | Pi |, and 
keep just the first term, which is of order \pi p. The resulting differential cross section is 

_ ^ )^/c~^ — Q,2^ -^2 ^ (23) 

do 2mc^ kp ^ 2mc^ c^kp 

We notice that this result is independent of the angle jS. For hQ = 10 eV, it is 24 orders of 
magnitude larger than the bare vacuum process (Fig.|^), whose total integrated cross section 
in this regime is known to be 0] 

a = 0.03a V2(^)6 , hn < 2mc'^ , 

where = e'^/mc^ is the classical electron radius. 

The opposite limit k ^ kp is applicable in the range hckp <^ HQ <^ mc^. For k not too 

large, and the angle (3 not too small, or too close to vr, we expand the Lindhard function 

n°(pj, Ui) in (p^) in powers of 1/ | | and keep just the first term, which is of order 1/ \ Pi p. 

For the differential cross section we then obtain 

da^ _ ^ ^ Ep kp Q 1 2 _ 4q;^ Ep 1 c^k% 

lb ~ 9^" 2^^~k' ^ ~ 97r2 2mc2sin^/3 • i 
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Fig.0 shows a log-log plot of ^ versus photon frequency Vl for /3 = 7r/2, together with 
the asymptotic dependences (p3D and (p^ . We see that the differential cross section has a 
maximum for frequencies corresponding to the photon wavevector k kp, where, with the 
present parameters, it is of order 10~^^ cm^. This is still a very small value, in fact somewhat 
smaller than the maximal photon-photon scattering cross section in vacuum. However, for 
photon frequencies in the ultraviolet and soft X-ray range, up to several keV's, this cross 
section is still about 15 orders of magnitude larger than that of the bare vacuum process 
(Fig.^). This is what may make it easier to measure in the end. 

So far we have dealt with the 2nd order scattering processes only. The 3rd order scattering 
amplitude ([14D gives rise to "triangle" diagrams, like that on Fig.|^. As we see, in these 
diagrams one also encounters the plasmon series, since one of the vertices corresponds to 
the interaction term Hi2. We expect the contribution of these diagrams to lie somewhere 
between that of the 2nd order processes, which we have investigated, and the 4th order ones, 
which we will discuss next. We shall not evaluate any of these 3rd order processes explicitly. 

IV. HIGHER ORDER PHOTON PHOTON SCATTERING PROCESSES: THE 
FOURTH ORDER (SQUARE) DIAGRAMS 

The straightforward second-order calculation of the previous section contains the main 
result of this paper, for the photon-photon elastic cross section in an electron gas. However, 
it is found, that for ^ ckp, the cross section has a very fast falloff with frequency, like 
much faster than the VL~^ behaviour known for fourth-order processes In fact, it 
can be anticipated, that higher order processes should become increasingly important as 
frequency increases. Therefore, it is at least in principle possible that third or fourth-order 
processes could become important when Vl is both very large with respect to ckp-, while 
still very small relative to mc^ (the optimal photon frequency for vacuum electron-positron 
processes). For this reason, we have undertaken the much heavier task of calculating the 
fourth-order scattering processes for two photons in the EG. A second reason for doing 
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this, is that this calculation will permit a better formal comparison with photon-photon 
scattering in vacuum, which is strictly fourth-order. 

We now consider the fourth-order scattering amplitude (|15|), which is the non-relativistic 
analog of the vacuum processes with the vertex corresponding now to the interaction term 
Hii- 

Substituting for Hn from ( p!o|) we get 



J dxt{%n\ a^^,agT[{A{x^).Vi) ... (l(x4).V4)] 4^4^ 
X I T[^^{x,)^{x,) . . . ^^{x,)^ix,)] I , (25) 

where each of the gradients Vj acts only on "^{xi). Evaluating this expression in the usual 
way, using Wick's theorem, we arrive at the already mentioned six diagrams (Fig.^), corre- 
sponding to the six possibilities for contraction of the operators. 
Making use of 

V,G\x, y) = (277)-^ j d'p me'^^^~y^G\p) , (26) 

and of relations (|T^) and (P^), and collecting together all the numerical factors we get finally, 
omitting all details, an expression for Mf 



diagrams 

We have denoted as the integral corresponding to the i-th diagram 

I, = J d^p (eli.pa) . . . {e,,.p,,)G%pa)G%p,2)G\p,,)G%p,,) , (28) 

where e*j„ is the photon polarization unit vector in z-th diagram, vertex n. 

Now since each electron Green's function G^ consists of two parts corresponding to 
electrons and holes, respectively, each of the above integrals splits into 16 terms. Each of 
these will contain a product of 4 step functions 0(±(| pj \ —kp)), which defines an integration 
region that is an intersection of interiors and/or exteriors of four mutually displaced Fermi 
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spheres. For general direction and magnitudes of photon wavevectors, the integration over 
such region would be quite intricate. For the sake of simplicity, we have not attempted 
to evaluate all integrals for a general situation. In the Appendix we present details of a 
particular calculation for the special geometry of forward scattering of two photons with 
opposite wavevectors (i.e. two initial photons with wavevectors k, —k scatter into the same 
two final photons), where the integrations turn out to be particularly easy. We shall just 
quote here the final results of this rather tedious calculation, where approximations have 
also been made based on relations (@) and (H). 

We have considered two limiting cases for the photon wavevector. First we take the 
case k ^ kp, which applies very well for ultraviolet photon energies HQ <^ 1 keV. We 



have obtained the result, Eq. ([67|), which, after passing from the rescaled variables (used 
throughout the appendix) back to the original ones, reads 



diagrams ^ 



where k and Q are the wavevector and the frequency of the photon. The function 
02, 01, ^2) is a geometrical factor, which accounts for the polarization of the photons, 
and is defined in equation (BSI) in the Appendix. 



Substituting for Q the bare photon frequency Q = ck, we get for M^f 

The corresponding fourth-order differential scattering cross section ^ is then finally, ac- 
cording to ([TSp 

To get an order of magnitude for ^ we substitute in the last expression the values typical 
for metals, i.e. kp ~ lO^cm"^ and Ep ~ 5 eV. With hVL ~ Ep, we obtain 

^ ~ (4 X 10-2) . (3 X 10-9) .(10-^)=^ . 1 . (10S)-2 ^ lO-^^cm^ . (32) 
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The second limiting case we have investigated is the case of photon wavevectors falhng 
above 2kp, which means X-ray photons with energies HQ >~ 1 keV. The sum over the 
diagrams is given by ([70|), which translated back to true variables reads 

diagrams ^ 

Substituting again the bare photon frequency for VL we obtain the scattering amplitude 

Mf = ^-^a\^^f^n-'E,k-,\ (34) 

and we notice that this does not depend on Vt. The differential scattering cross section then 
equals 

We estimate the order of magnitude of ^ for k 2kF and obtain 

^ ~ (4 X 10-2) . (3 X 10-9) . (10-^)^ . {10~'^y . (10«cm-i)-2 ^ IQ-^^cm^ . (36) 
do 



Equations ( pi]) and (^) represent the main results of this section. We can now compare 
them with their equivalent results in vacuum, namely ||^ 

a = 0.03aV^(^)6, nVL<^2mc^ (37) 
cr = 4.7a^(^)2, nVL^2m^, (38) 

where a is the total integrated cross section for unpolarized photons and rg = e^/mc^ is the 
classical electron radius. We note, first of all, that the dependence at low frequencies is 
the same as in our low frequency second-order result of Section 3, Eq. (p3D , while the high 
frequency behaviour is the same as in both the second-order result of Eq. (p^) , and in 
the present fourth-order results (|3T| ) and (pS]). Secondly, we see that the Qr^ frequency 
dependence (assuming fii = ^2 = ^) is the same for both 4th order in vacuum and 4th 
order in the EG, provided the high-frequency regime is reached in each case. Numerically, 
however, the coefficient in front of $7"^ is much smaller in the EG case. At lower frequencies, 
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instead, we find a new regime in the EG case, where the 4th order cross section falls off like 
Eq.(^). This regime does not exist in vacuum, and is clearly due to the presence of a 
Fermi surface in the electron gas problem. 

We can now try an overall graphical comparison of all these results. This is sketched in 
Fig.^. We see that the 2nd order processes in the EG are dominant for photon energies in the 
range from uop up to ~ 10^ -t- 10^ eV. For X-ray photons of 10^ eV the relativistic vacuum 
processes become important and dominate at all higher photon frequencies. Hence, the 4th 
order EG contributions are always negligible, and by a huge factor, in spite of their slower 
falloff with photon energy. It is reasonable, therefore, to assume without proof that similar 
conclusions will apply to the third order processes, which we have accordingly ignored. 

V. DISCUSSION, AND POSSIBLE EXPERIMENTAL VERIFICATION. 

We have considered the problem of elastic ultraviolet and X-ray photon-photon scattering 
inside a solid, idealized as a free electron gas. We find that the presence of the electron gas 
should give rise to new scattering processes, much more important than those present in 
vacuum. There, only 4th order processes, important for 7-ray photons in the MeV range, are 
operative. In the nonrelativistic electron gas, instead, the existence of two-photon vertices 
introduces large second order processes, which are most efficient for photon wavevector 
roughly equal to the Fermi wavevector of the metal. We also predict important plasmon 
resonances to take place, and an angular dependence different from that of the vacuum 
processes. 

We have also recalculated the 4th order processes in the EG. We find that these are 
themselves different from those in vacuum, mostly due to the breakdown of electron-hole 
symmetry, and their importance increases as the photon frequency decreases. However, the 
plasma frequency up represents in practice a lower bound for the frequency of a photon, if it 
should penetrate inside a metal, and we have found the 4th order processes to be negligible 
for frequencies higher than ujp. The overall situation is summarized by Fig.|^, which also 
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gives an order of magnitude of cross sections for a metal such as potassium, chosen because 
of its low plasma frequency. 

Let us now briefly consider the possibilities for experimental verification of our calculated 
electron-mediated photon-photon elastic scattering processes. We shall discuss two cases, 
both extremely idealized. In the first we consider the use of a powerful pulsed laser source, 
operating in the ultraviolet region, hQ ~ 10 eV. The second case will be that of an X-ray in 1 
keV range, such as that, which can be obtained from a synchrotron source (plus undulators). 

In order to estimate the scattering cross section for the first case, we refer to Fig.^. 
Making use of a plasmon resonance, we can expect ^ ~ 10"^^ cm^ for a metal. One could 
think, for example, of a crossed-beam experiment. Two laser beams are focused on a small 
area and penetrate into the metal and the scattered radiation is picked up by a detector. 
Without worrying about problems of power dissipation, let us think of one laser pulse per 
second, of a duration of T picoseconds and a peak power of Jx 10^° Watts. The pulse contains 
an energy of ~ I T x 10^^ J, which corresponds to ~ I T x 10^^ photons with energy 10 eV. 
The spatial length of the pulse is 3 T x 10^^ cm = 3 T x 10^ A. Imagining a penetration 
depth of ~ 200 A, the number of photons inside the metal is / x 10^^. The total scattering 
cross section, as seen by another photon, is therefore ~ / x 10"^'^ cm^. If the second beam, 
having the same parameters as the first one, is focused on the area of (3 /im)^ ~ 10~^ cm^, 
the number of collisions during the pulse is T X I X T X 10^^ ~ P T. For 

a peak power of, e.g., 10^° W (/ = 1) and a pulse duration of 1 ps (T = 1) we can expect 
about one scattered photon per unit solid angle per second. This flux should be, presumably, 
detectable with present experimental techniques. 

In the case of a synchrotron source, one can work in the soft X-ray region, where the 2nd 
order scattering cross section has a maximum. According to Fig.|^, this occurs for photon 
energies h^l ~ hckp in the keV range, where one has ^ ~ 10^^2 

cm . The penetration 

depth in the metal in this region might be ~ 10^ A. To be more specific, we consider a 
beam of 1 keV radiation, consisting of 10^ pulses per second, with pulse duration 100 ps and 
peak photon flux 10^° photon s~^. If the cross section of the beam is taken ~ 10~^ cm^, the 
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resulting number of scattered photons is again ~ 1 photon per second and unit sohd angle, 
which is similar to the crossed UV laser case, and should have no sample burning problems. 

In all cases, it is clear that, apart from the products of photon-photon collisions (which 
we are after) there will be a background of photons scattered by other entities due to different 
mechanisms. There will be, for example, single-photon processes, like photon scattering from 
free electrons. The order of magnitude of the corresponding cross section can be estimated 
from classical Thomson formula to be of order (ia ~ ~ 10~^^cm^, and is in general larger 
than the photon-photon cross section, at least for reasonable intensities. These single-photon 
processes, however, should be distinguishable from the two-photon processes of our interest, 
either by requiring photon-photon coincidence, or by making use of either the quadratic 
intensity dependence of the photon-photon events, or of their polarization dependence (p2D. 

A subtler complication is that of additional contribution of interband transitions to two- 
photon processes. Since interband transitions are vertical, the intermediate states in the 
scattering process can be real ones, at particular photon frequencies. Care should be taken, 
therefore, to stay away from these frequencies as much as possible. As for inelastic processes, 
which probably represent a big problem, one may perhaps think of using the typical elastic 
polarization dependence obtained in Eq. (P^), to subtract them out. 

Yet another possible route, although admittedly very speculative, could be studying the 
subtle changes in the interference pattern which one can expect to take place at very high 
photon intensity in a classical interference experiment. The classical theory of interference 
is of course based on the superposition principle, i.e. on photon-photon scattering being 
exactly zero. If that is no longer true, one can expect that correlation effects building up 
between the photons should in principle change intensities in the interference pattern. De 
Martini et al. |T0[ , for example, have recently shown how random selection of polarization in 
front of the interferometer presents precisely a realization of this kind of effect. In that case, 
when the beams are very intense, photons become correlated with one another through the 
random polarizer. As a result, their distribution tends to become Bose-like, which can be 
imagined as "trains" of photons going down one or the other slit separately. This has been 
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shown to weaken the interference pattern in a characteristic and measurable manner. 



In principle, when the beam is sufficiently intense and the "slits" are metallic (e.g., a half- 
silvered mirror), a similar kind of modification of the interference pattern should be expected 
even without the random polarizer. Our estimates of the cross section could in principle be 
used to evaluate the threshold value of the intensity in a well-defined interference geometry. 

A general remark may be in order, before closing this paper. As things stand, there are 
two elements which suggest that our calculation might remain, at least for while, of purely 
academic interest. The ffist element is that our calculated photon-photon cross sections are 
still very small. Detection of these processes may require a nontrivial experimental effort. 
The second point is that detection of our proposed process does not in itself provide new 
basic information, on either the photon, or the system. 

Our viewpoint on these issues is open. Without embarking in a discussion of why one 
should or should not try to measure the process we calculated, we have simply meant to 
point out, for future record, that there are new photon-photon processes, which were not 
discussed so far, and which are in principle of measurable intensity. 

APPENDIX: DETAILS OF FOURTH ORDER CALCULATION. 

To avoid problems with ill-defined expressions we shall evaluate Mjf for the case of 
scattering at a small angle, i.e. ki = k, k2 = —k, k[ = k + q, k'2 = —k — q, where q must be 
such that I k + q\ = \ k |, and then perform the limit | g |— to obtain the forward scattering 
amplitude. Because of the transversality of the vector potential, in the limit | g |— > we 
have 

k.ei = k.e2 = k.e^ = fc.Cj = . (39) 

In order to deal with dimensionless variables we rescale quantities in the integrals (|28|) and 
introduce dimensionless frequencies and wavevectors by the relations 

.^M., (40) 
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k = kpk' . (41) 

The integrals l'^ calculated with the rescaled variables lu' , k' are related to the original ones 
by 

h = (^fkFl[ , (42) 

and in the following the primes in all dimensionless quantities will be omitted. 

As already mentioned, each of the integrals (p8|) splits into 16 terms corresponding to 
all possible combinations of electrons and holes participating in the process. Some of these 
terms are, however, immediately seen to be equal to zero. Apart from obvious vanishing of 
the terms in which all 4 internal lines correspond simultaneously to electrons, resp. holes, 
terms in diagrams 1 and 3 in which there is an electron and a hole in the state with the same 
3-momentum, namely p2 = p^-, have to vanish, too, due to the vanishing of the integration 
region. In the Tab. | we enumerate all possibly non-zero contributions to integrals Ji, h-, I3, 
corresponding respectively to the diagrams 1,2,3 (Fig. p!o|) , where + sign denotes a hole and 
— sign an electron. 

Each of the terms has a general form 

/d.K,so ... ... 5!±%^. (43) 



where the momenta pn, . . . ,Pi4 depend on p and frequencies Un, . . . , depend on u. Per- 
forming the frequency integrals using the residue theorem and passing to the limit | ^ | 
wherever this can be done in a straightforward way we get the following expressions 

ha + hb = ^7ri J d%..]&{l-\p\)e{\p-k\-l) 

-Ak.p + 2k'^ 

X =; (44) 

{n + 2k.p-k^y{n-2k.p + k^y 

ha + hb = 47rz I d% . .] 9(1- I p 1)9(1- \p + k\)Qi\p-k\ -1) 

X 4 i (45) 

Ak.p {n - 2k.p + k^y 

ha + hb = ^rcz J d%..]Q{l-\p\)Q{l-\p + k\)Q{\p-k\-l) 
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1 



1 



1 



X 



Ak.p n - 2k.p + F -n - 2k.p + P 

/2d + /2e = 47rz J d%..]e{\p\ -i)e{i- \p + k\)e{\p~k 
1 1 



X 



Ak.p {n + 2k.p + k'^y 

hd + he = 47rz J d% . .] 0(1 p I -1)6(1- I p + k 1)6(1 p- k 
1 1 1 



-1^ 



-1) 



X 



Ak.p n + 2k.p + P -fi + 2k.p + F 



(46) 



(47) 



(48) 



he = -27rz / d% . .] 6(1- I p 1)6(1 p + k I -1)6(1 p- k \ -1) 



X 



4k.p 



{Q + 2k.p + k^y {Q-2k.p + k^y 



he = -2111 / d% . .] 6(1- \p\)ei\p + k\ -1)6(1 p- k I -1) 



X 



2Q 



1 



1 



1 



1 



(49) 



(50) 



n + 2k.p + k^ n- 2k.p + P n- 2k.p- k"^ VL + 2k.p- P. 
where VL is the photon frequency. We have used the relations (|39|) and introduced the 
notation 



(eii-p) {ei2.p) {e,s.p) {e^.p) = (ei.p) {e2.p) (e^.p) (e^.p) 



(51) 



in the above expressions. In several terms we have performed a shift and/or an inversion of 
the integration variable p. Terms 2f and 3f are identically zero since the corresponding 6 
functions have zero product. 

Terms he - hf and hg - hn require a comment. They contain products like 6(1— | p \ 
)6(| p + q\ —1), which restrict the integration volume to a thin shell about a half of the 
surface of the Fermi sphere (see Fig.|TT]). The volume element of the integration region can 
then be written as 

dV = ±dS.q = ±dSp.q. 
At the same time these terms contain one denominator of the form 



{p + q) —p =2p.q + q ^ 2p.q 
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which in the hmit | ^ precisely cancels the same product in the volume element. In 
the limit | > therefore these terms become well defined surface integrals. Summing all 
the contributions from the diagrams 1, resp. 2, we obtain expressions 



hs — he + hd + he + h 



f 



= -TTi[ dS[...] i 1 (52) 

Ap\=i n + 2k.p- -n + 2k.p- 

hs = hh + hi + hj + hi + hm + hn 

= ~mf dS[...] i i , (53) 

^|p1=i Tl + 2k.p + k^ n-2k.p + k^ 

where the integration region is the whole surface of the Fermi sphere. Terms hg and hk are 

zero again because of vanishing of the corresponding product of functions. 

Now we must take into account also contributions of other 3 diagrams which differ from 



those on Fig.|TO| by the orientation of the internal loop. Due to the absence of charge- 
conjugation symmetry, Furry's theorem does not apply, and the contributions of diagrams 
which differ from each other by the orientation of the internal loop may not be equal. It can 
be easily seen that the contributions of diagrams 1 and 3 do not depend on the orientation, 
while in case of diagram 2 the change of the orientation is equivalent to the change of sign 
of the photon frequency Q. Therefore the sum of contributions of all the diagrams can be 
written as 

J2 h = ^{ha + hb) + 2/i, + 2{ha + hb) + 2/3C + 2{hd + he) 
all diagrams 

+ ha{^) + ha{-^) + hb^ + hb{-^) + he^ + /2c(-^^) 

+ hdm + hd{-^) + /2e(^^) + he{-^) + hsm + hs{-^) • (54) 

At this point we shall perform some approximations. First we make use of the conditions 
(0), d^), which we assumed at the very beginning of our non-relativistic treatment. We write 
the bare photon dispersion relation Vt = ck m the rescaled variables as 

VL = Ak, (55) 

where 
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Substituting to this expression a typical Fermi wavevector of metals kp ~ 10*^ cm~^ we 
observe that 

A ~ 10^ > 1 (57) 

holds for A. Actually, the last condition is just the square root of the condition (^. More- 
over, the condition (|^) now turns out to be equivalent to 

k<t:A. (58) 

Because there is always a hole participating in the process, the integration variable p 
must fall inside at least one of the 3 Fermi spheres defined by|p|<l, |p + /c|<l and 
\ p — k I < 1 . This implies 

\p\<l + k, (59) 

and due to (0) and (|58D also 

\p\<^A. (60) 

Now in the denominators of our integrals we encounter terms like Q ± 2k.p ± k^ = 
Ak ± 2k.p ± fc^, and the relations (|58D, ( |60D and ( [57| ) immediately suggest that we can 
expand the integrands in powers of 1/A and keep just the first non-zero term. Since the 
expression (|5^) is an even function of Q and therefore of A, the expansion contains only even 
powers, and represents actually an expansion in powers of ^ = The first non-zero 

term turns out to be of the order 1/A'^ in all the integrands, and neglecting terms of the 
order of 1/A^ we obtain respectively 

2{ha + h,) = TTZ I ciM- • •]0(1- I P 1)0(1 P-k I (61) 
2(/3a + hh) + {ha + hh)m + {ha + hh){-^) = 
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« J d% . .] 6(1- I pmi- I p + k 1)6(1 p- k I -l)^L_P__± (62) 

2(/3d + /Se) + {hd + he^) + {hd + /2e)(-^^) = 

vr^ j d'^. . .] 6(1 p I -1)6(1- I p + k 1)6(1 p- k I _i)i_ (2fc.p + P)^ ^^^^ 



k.p 

ic = 2h, + i2cin) + hci-n) 

r ^ ^ —}6k^ 

= mj d% . .] 6(1- I p 1)6(1 p + k I -1)6(1 p - I -1)-^ , (64) 

is = 2h, + hsin) + i2si-n) 

= m dS[..]——{2k.p + k^), (65) 

j\ p\=i \r 

where we have aheady grouped together some terms. 

Now we could, in principle, calculate the above integrals in a straightforward way and 



obtain the scattering amplitude for any fc, restricted only by the condition ([58|). For sake of 
simplicity we shall do the integrations explicitly for 2 limiting cases, in which the integration 
region becomes particularly simple, first for A; <^ 1, and then for 2 < A; ^ A. 

The condition A; <^ 1 allows us once again to pass in several terms to surface integration 
(in the same way as above when we were performing the limit | g* |— 0). In the limit A; <^ 1, 
the integration regions of the terms la,lb,2a,2b,3a,3b,2d,2e,3d and 3e become the same and 
equal to the half of the surface of the Fermi sphere 

I p 1= 1, for which k.p < holds. Collecting these terms together and denoting their sum 
as // we obtain 

// = 2{ha + hh) + 2{ha + hb) + {ha + hhW) + {ha + hb){-^) 
+ 2{hd + he) + {hd + /2e)(^^) + {hd + he){-^) 

= ^ L dS[. . .] {8k^ + IQPk.p) . (66) 

il Jk.p<Q 

The remaining terms in (0) are the volume integral Ic (^^^ ^^"^ surface integral Is (§^ . 

Since we are considering the limiting case A; <^ 1, it is now enough to identify in the sum 
h + Ic + Is terms which contain the leading power of k. In // the leading term in k is of 
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order k^. It is easily seen that the integration volume of the integral Ic behaves as and 
the contribution of this term is therefore of order k^. In the surface integral Is we encounter 
a term ~ k'^k.p, which, however, vanishes after the integration over the whole surface of 
the Fermi sphere and the remaining term contributes as k^. In the limit ^ 1 therefore 
the dominant contribution to the scattering amplitude comes from the leading term in (pBD, 
which is readily evaluated as 

E /. = ^^Wi,02,0u;). (67) 

all diagrams 

In the last expression we introduced a function 

/ / /'27r ^ / TT 

5(01,02,01,02)= / (i0COs(0 - 0i) COS(0 - 02) COS(0 - 0,) COS(0 - 02) = - X 

Jo 4 

X COS(01 + 02 - 01 - 02) + COS(01 - 02 + 01 - 02) + COs(0i - 02 - 0i + 02) , (68) 

which represents a factor taking into account the polarization of the photons, and 
01, 02, 01, 02 cire direction angles of the unit polarization vectors of the photons. 

In another limiting case, 2 < k A, the Fermi spheres are displaced from one another 
by a distance larger than 2 radii, and therefore do not intersect anymore. The integration 
region of integrals l2a + and I^a + -^36 becomes empty and that of the remaining volume 
integrals becomes the whole interior of one of the spheres. Denoting the sum of all volume 
integrals as ly we obtain 



Iv = ^Jd%..]Qil-\p\) 



^ 4(2A;.p-P)2 
k.{p — k) 



X 5(01, 02, 0'i, 0^) {^(-33 + - 15k^) + he - Se + ?,e - e) in • (69) 

45 o fc + 1 



Apart from this, there is also the surface integral Is (|65D, which equals 

/. = -^^i?(0i,02,0'i,0;). (70) 

Now for large k, the leading power of k in Jy is k'^. The dominant contribution to the 
scattering amplitude will therefore come from Is, which behaves as k'^. However, on the 
side of large k we are still restricted by the condition (pSf ), and the validity of expressions 
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(|69|) and (^) is thus limited, say, to 2 < < 10. For this interval of /c, we have numerically 
checked the relative contributions of ly and Is, and really found Jy to be negligible. The 
final result is therefore given by (|70[). 
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FIGURES 



FIG. 1. The vacuum photon-photon fourth-order scattering diagrams. 
FIG. 2. Second order photon-photon scattering processes in a nonrelativistic electron gas. 

FIG. 3. A diagram with zero ampUtude. 

FIG. 4. Screening of the bare polarizabihty in the second order diagram. 

FIG. 5. Photon-photon planar scattering geometry (chosen for simplicity of illustration), ki, ^2 
incoming photons, ki,k'2 outgoing photons. 

FIG. 6. Differential cross section versus angle /3 for the 2nd order processes in the planar 
geometry of Fig. ^, and for 6 = tt/2. The peaks are due to the plasmon resonance. Their 
sharpness is a consequence of the random-phase approximation, where the plasmon lifetime is 
infinite. 

FIG. 7. A log-log plot of the calculated differential cross section versus photon frequency for 
the 2nd order elastic photon-photon scattering processes in an electron gas (away from plasmon 
resonances). Parameters chosen are representative for potassium metal. 

FIG. 8. "Triangle" diagrams corresponding to the 3rd order processes. 

FIG. 9. An overall graphical comparison of the photon-photon elastic cross section in an elec- 
tron gas and in vacuum. In an actual solid, or liquid, both contributions will be present. 

FIG. 10. Three diagrams corresponding to the 4th order processes in an electron gas for the 
particular case of the scattering of two incident photons with opposite wavevectors. 

FIG. 11. The integration region, which is a thin outer shell of the Fermi sphere. 
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TABLES 



Diagram 1 Diagram 2 Diagram 3 
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TABLE I. Terms giving a non-zero contribution to the diagrams 1-3. + denotes a hole, — an 
electron. 
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